Relaxation energies and excited state structures of poly(para-phenylene) 
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We investigate the relaxation energies and excited state geometries of the light emitting poly- 
mer, poly (para- phenylene). We solve the Pariser-Parr-Pople-Peierls model using the density matrix 
renormalization group method. We find that the lattice relaxation of the dipole-active l 1 B± u state 
is quite different from that of the i S B^ u state and the dipole-inactive 2 1 A~^ state. In particular, 
\Q . the l 1 ^!^ state is rather weakly coupled to the lattice and has a rather small relaxation energy 

£Nj ■ ~ 0.1 eV. In contrast, the 1 3 B^ U and 2 1 A£ states are strongly coupled with relaxation energies of 

~ 0.5 and ~ 1.0 eV, respectively. By analogy to linear polyenes, we argue that this difference can 
be understood by the different kind of solitons present in the I 1 B^ u , 1 3 -B^„ and 2 1 states. The 
D ' difference in relaxation energies of the l 1 B± u and l 3 B^ u states accounts for approximately one-third 

t 1 , of the exchange gap in light-emitting polymers. 

& ■ 



PACS numbers: 71.10.Fd, 71.20.Rv, 71.35.Aa 



I. INTRODUCTION 

T3 

Electron-lattice coupling has profound effects on the behavior of conjugated polymers. It is responsible for the self- 
trapping of excited states and plays a vital role in determining the interconversion between excited states. Predicting 
interconversion rates is important for understanding many electronic processes in conjugated polymers, e.g. the 
determination of the singlet exciton yield in light emitting polymers. Although not (directly) addressed here, the 
formation of quinoid-like structures in poly(para-phenylene) is largely responsible for the tendency of these molecules 
to planarize in the excited state. 

Electron-electron interaction is also important in conjugated polymers. It is responsible for excitons with large 
binding energies and a large singlet-triplet exchange gap. The inter-play of both electron-lattice and electron-electron 
interactions leads to a complicated, but interesting description of the excited states. In this paper we use the density 
t-H ' matrix renormalization group (DMRG)* method to solve the Pariser-Parr-Pople-Peierls model - a 7r-electron model 
that treats both electron-lattice and electron-electron interactions. 

Poly(para-phenylene) has a l x A^ ground state. The lowest lying optically allowed state is the l 1 !!^ state, and 
since this is the emitting state in LEDs, its behavior is obviously important. Another important state is the I 1 A+ 
state. In light emitting polymers it is argued that this is the rn^Ag stated, which is strongly dipole connected to the 
1 -Bi„ state and consequently plays an important role in non-linear spectroscopy. As a consequence of both strong 
electronic correlations and electron-lattice coupling in linear polyenes, the relaxed energy of the 2 1 A+ state lies below 
the relaxed energy of the l x B^ u state, rendering these systems non-electroluminescent. We investigate the behavior 
of this state in poly(para-phenylene). Because charge carriers are injected with arbitrary spin, triplet formation is 
also of importance, so we also explore the behavior of the 1 3 B^ U state. 

The phenyl-based conjugated polymers are extrinsically semiconducting as a consequence of the chemical structure 
determined by the a bonds in the absence of 7r-conjugation. Thus, with all bond lengths equal there is still a 
semiconducting band gap. However, as for linear polyenes with extrinsic dimerization, coupling of the 7r-electrons 
to the lattice is still important, as it causes bond lengths to change resulting in exotic types of bound non-linear 
excitations. Understanding how the non-interacting description of the excited state structures change as a consequence 
of electronic interactions is a key goal of this work. To facilitate this goal, we map the excited state structures of 
poly(para-phenylene) into the more familiar description of mid-gap single-particle states and the associated solitons 
of linear polyenes. 

Electron-lattice coupling in light emitting polymers has been investigated by a number of groups using a variety 
of methods. Beljonne et al. studied the relaxation of the 1B U singlet and triplet in short poly(paraphenylene 
vinylene) oligomers by solving intermediate neglect of differential overlap models using multi-reference configuration 
interactions. Although they found larger relaxation energies for the triplet than the singlet, the difference of ca. 0.1 
eV in the four-ring oligomer is considerably smaller than our results as presented in section IV. B. Ambrosch-Draxl et 
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alX performed a density functional theory calculation on the ground state neutral structure of PPP, while Zojer et 
al. 5 performed semi-empirical Austin Model 1 calculations on the ground state neutral and charged structures of PPP. 
Finally, Artacho et al& performed a GWA-Bethe-Salpete equation calculation on the state of PPP. We compare 

their predictions to ours in the results section. To our knowledge, this work presents the first large-scale calculation 
of the relaxation energies and geometrical structures of the l 1 i3 1 ^ u , 1 3 B^ U and 2 X A+ states in poly(para-phenylene) 
oligomers. 

The experimental relaxation energies of the lowest-lying singlet exciton and doped polaron in a wide variety of 
polymers of different conjugation lengths has recently been presented by Wohlgenannt 7 . Other experimental results 
are discussed in the results section, IV. B. 

Extensive investigations of the excited states of poly(para-phenylene) for fixed geometries within the Pariser-Parr- 
Pople model, solved by the DMRG method, are described in re&. 

In the next section we introduce the Pariser-Parr-Pople-Peierls model of conjugated polymers. Next, we solve 
the Pariser-Parr-Pople-Peierls model in the non- interacting limit, namely the Peierls model. We discuss soliton 
wavefunctions and soliton-antisoliton confinement. Then we solve the full interacting model with the DMRG technique. 
We briefly discuss the DMRG algorithm before describing the results. Finally, we summarize and conclude. 



II. THE PARISER-PARR-POPLE-PEIERLS MODEL 



The Pariser-Parr-Pople-Peierls model, Hpppp, is a tight-binding model of the 7r-electrons that includes both long 
range Coulomb interactions and electron-lattice coupling. The electrons and lattice are coupled together by the effects 
of changes in the bond lengths both on the one-electron transfer integrals and the Coulomb interactions. We treat 
these effects up to first order in the change of bond length. As the density-density correlator, (Ni — l)(Nj — 1), decays 
rapidly with distance^, it is also a reasonable approximation to retain changes in the Coulomb potential for only 
nearest neighbor interactions. 

We thus define the Pariser-Parr-Pople-Peierls model as, 

Hppp P = -2j2ufi + Wj2MNi + i-l)(Ni-l) (1) 

i i 



f < = J E (^1,^+4^+1,.) (2) 



is the bond-order operator, where c\ a creates an electron with spin a in the 7r-orbital on site i. The one-electron 
transfer integral is 

U = t + ^- 1 (3) 

where 

A 4 = -2a(u i+ i - Ui), (4) 
and Ui is the displacement of the ith atom from its equilibrium position. Thus, we define, 

r _ SU A, 

Su t = = - — , (5) 

a 2a 

as the change in length of the ith bond from its initial - undistorted - value (determined by the a bonds) . 
The electron-phonon coupling constant is, 
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and Vij is the Ohno potential, 

U 

ij ~ A /l + (t^-/14.397p' (7) 

for the reference (that is, the undistorted) structure, where the bond lengths are in A. We set U — 10.06 eV, t = 2.514 
eV, A = 0.12, a — 4.67 eVA~ 2 and the undistorted bond length as 1.405 A. The parameters used are essentially those 
of reS, with t and a chosen such that if the bond lengths are those of reS, the transfer integrals are as well, while A 
was chosen to obtain good agreement with the ground state structure used in refi. 

The second term on the right-hand-side of Eq. Q is the change in the Coulomb interactions from changes in bond 
length, where 

Wd^\ Ur (U/U.397r 

2a \dTmn) Tti=ro 2a(l + (C/r /14.397) 2 )3/2- y > 
It is instructive to re-write this term as 

- 2aWj2(^+i - u t ){N l+1 - 1)(JV, - 1), (9) 

i 

where we have used Eq. @. Expanding and resuming we see that this term has two components. One component is 
the electron-phonon coupling arising from the change in the ionic potentials, 

-2aW^{u i+1 -u i - 1 )N i . (10) 

i 

The other component represents the changes in the nearest neighbor electron-electron interaction from the change in 
bond length, 

- 2aWj2(u i+ i ~ Ui)N i+1 Ni. (11) 

i 

The first term on the right-hand-side of Eq. is just the electron-phonon coupling arising from the change in the 
kinetic energy. 

Notice that this model does not describe free rotations of phenyl rings relative to one another. Thus, its applicability 
is to ladder poly(para-phenylene), where the stereo-chemistry causes the rings to have a planar geometry, or polymers 
in the solid state, where ring rotations are more restricted. 

As described in previous paper a 9 i 10 , by using the Hellmann-Feynman theorem we can derive a self-consistent 
equation for A, for any state, 

A; = 27rfA (<T;) - W{Di) - r) , (12) 

where 

Di = (N i+1 - l)(Ni - 1) (13) 

is the density-density correlator for the ith bond. T is determined by imposing constant chain lengths, Aj = 0, 
implying that, 



T={T i )-W{D i ), (14) 
where the over-bar represents the spatial average. 



III. NON-INTERACTING LIMIT 



In this section we describe the solutions of the Pariser-Parr-Pople-Peierls model in the non-interacting limit (U = 0), 
namely the Peierls model. 

Fig- IH shows the fractional change in transfer integrals for the ground state, 5ti, defined in Eq. JSJ). The bonds are 
defined by Fig- EI Since the bonds are initially all of the same length, we see that the coupling of the 7r-electrons to the 
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FIG. 1: The fractional change in transfer integrals of poly(para-phenylene) from the uniform value, t, in the non-interacting 
limit. The relevant electronic states are indicated in the panels. The labels refer to the bonds shown in Fig. [5] Only the upper 
rung of bonds are shown. Notice that the change in transfer integrals is opposite to the change in bond lengths. 



lattice has caused an effective 'bond' alternation. The phenyl-ring bonds shorten while the bridging bond lengthens. 
This is the benzoid structure, as the phenyl-ring bonds are roughly all of the same length. 
To see this effective bond alternation we define the summed bond distortions as, 

St n = Sti; odd n (15) 

iephenyl ring 

and 

5t„ = 5t i= bridging bond ; even n. (16) 
Then we defined the normalized, staggered and summed 'bond' alternation, 5 n , as, 

*n = ^(-l) n . (17) 

Fig.yjshows S n for the ground state. Under this mapping the phenyl ring is equivalent to a double bond (or dimer) 
and the bridging bond is a single bond. As for polyenes, this effective alternation increases the semiconducting band 
gap. Note that end-effects coupled to the constraint of an overall constant contour length causes the oscillations in 5 n : 
there are greater distortions in the phenyl rings at the end of the chain than in those in the middle of the chain. Thus, 
in the middle of the chain the summed distortion in bond lengths in a phenyl ring is not quite equal and opposite to 
the distortion of the bridging bonds. 

Next consider the lBi u excited state structure, shown in Fig. ^ This is the quinoid structure, illustrated in Fig. [21 
(b). In contrast to the ground state, there is now a significant variation in the bond lengths in the phenyl-ring: bonds 
labelled 1 shorten, while bonds labelled 2 lengthen. The bridging bond also shortens. 

At first sight the excited state lattice distortions of poly(para-phenylcnc) represented in Fig.^Jdo not resemble that 
of a linear polyene. We therefore might enquire whether the geometrical defects (for example, solitons, polarons, etc.) 
and their associated mid-gap electronic states also exist in an analogous manner in poly(para-phcnylcnc). To show 
that bond defects do exist in an analogous manner to linear polyenes we again consider the summed bond distortions, 
defined by Eq. (|T7|). 

The 1-Biu state structure is illustrated in this way in Fig. EI The relaxed 1B\ U state creates a 'polaronic' structure, 
whereby the average bond length in the phenyl ring increases while the bridging bond length decreases, but there is 
no reversal in bond distortions from the ground state. This polaronic structure of excited states occurs in extrinsically 
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FIG. 2: (a) The bonds illustrated in Fig.^and Fig. El an d bond lengths in A of the ground state determined in the interacting 
limit, (b) The quinoid structure of the l 1 B^ u state. Bond lengths in the center of the distortion in the interacting limit. (See 
section ITV Bl and Fig. |SJ 
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FIG. 3: The staggered, normalized and summed bond distortions of poly(para-phenylene) (as defined in Eq. JEJ) in the 
non-interacting limit. 



semiconducting polymers where the ground state is non-degenerate: reversing the sign of the bond distortions gives a 
higher energy^i. A bond defect, or soliton, separates two regions of opposite bond distortions. Creating a soliton and 
antisoliton pair and moving them apart creates a region of reversed bonds. Thus, there is a linear confining potential 
between the soliton and antisoliton for large separations. As in linear polyenes, these bond defects are also associated 
with mid-gap states. 

Associated with the two- mid gap single-particle states of the excited state are a bonding, ipf , and anti-bonding, ip~ , 
molecular orbital (where i is a site index). These molecular orbitals are analogous to the bonding and anti-bonding 
orbitals of molecular hydrogeni£. The molecular orbitals are constructed from localized Wannier functions, 4>i and 
<j)i, which represent the soliton and anti-soliton respectively. In particular, 

^ = -^±&), (18) 

or inverting, 



6 



0.3 



0.25 



0.2 



-s- 0.15 



0.1 



0.05 



■ □ ■ □ 

5 10 15 20 25 30 35 40 



Phenyl ring, n 



FIG. 4: The soliton (solid symbols) and anti-soliton (open symbols) probability densities (defined by Eq. ^H ). 



and 

h = ^Pt -i^i)- (20) 

In linear polyenes with degenerate ground states the soliton and antisoliton are widely separated. However, as de- 
scribed above, they are confined in extrinsic semiconductors. This confinement is illustrated for poly(para-phenylene) 
in Fig. 01 which shows the soliton and antisoliton probability density summed over each phenyl ring, 

€= E ( 21 ) 
iephenyl ring 

We see that the soliton and anti-soliton wavefunctions are centered on neighboring phenyl rings in the middle of the 
chain. 

We now use this description of the molecular orbital defect states (Eq. 1181) to describe the excited states. As first 
shown by Ball et al. for trans-polyacetylenei^, as a consequence of spin and spatial symmetries the relaxed l l B^[ u 
and 1 3 B^~ U states have quite different solitonic characteristics. These differences become important when the spin 
degeneracy is lifted by electronic interactions, and they help explain the quite different geometrical distortions of 
these two states in the interacting limit. 

We first review the argument of re£i£. First, let us consider the singlet, l 1 -Bi, 1 state. We write this as, 

\l 1 B lu ) = ±(cl,cl l -cl l c%)\V), (22) 

where \V) represents the occupied sea of valence states and Cj_ ff creates an electron with spin a in the mid gap state 
|V' ± )- If c\ and cj. creates an electron in the states \4) and respectively, then 



4 - J_(J + -V 
V2 



c± CT = — ($,±4). (23) 



Inserting Eq. into Eq. we have 

ll 1 ^) = -U4eJ -c\c\) | V). (24) 



V2 

;atively cha 

no electrons, so it is positively charged and also spinless. Similarly, ejej creates a pair of electrons in the anti-soliton, 



ejej creates a pair of electrons in the soliton, so it is negatively charged and spinless, while the anti-soliton contains 



7 



TABLE I: The relaxation energies of the 1Bi u , 2A g and charged (polaron) states for para-phenylene oligomers (in eV) calculated 
from the Peierls model. 



Number of phenyl ring; 


5 


2A g 


polaron 


4 


0.23 


0.14 


0.06 
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0.10 


0.04 


20 
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0.02 


40 


0.06 


0.03 


0.01 



while the soliton contains no electrons. The l^Biu state is therefore a linear superposition of spinless positively and 
negatively charged soliton-antisoliton pairs. 

A similar argument applies to the triplet, l 3 Bi u state, 

|l 3 B lu > =±(cl^_ l +cl l cl r )\V) 

= -3(4 S ; + 44)^>' (25) 

showing that it is a linear superposition of neutral spin- 1; soliton-antisoliton pairs. 

Electronic interactions have a significant affect on these states. The spinless oppositely charged solitons of the 
^.Biu state bind to form an exciton-polaron. Conversely, the neutral spin-1/2 solitons of the 1 3 B\ U state do not bind, 
but cause a locally strong lattice distortion. We investigate these structures in the next section. 

To aid in our understand of exciton-polaron structures (to be described below) and relaxation energies, we also 
investigate charged (polaron) states in the non-interacting limit. Fig.|3shows the polaronic structure associated with 
a doped particle. Table [I] lists the relaxation energies of the lB lu , 2A g and charged (polaron) states for different 
oligomer lengths. We note that the relaxation energy of the 1Bi u state is considerably greater than for the polaron, 
and that the relaxation energies reduce as the oligomer lengths increase. 

In general, poly(para-phcnylcnc) is not planar because of the steric repulsion of the hydrogen atoms on neighboring 
phenyl rings. The torsional angle between adjacent phenyl rings for a single chain is estimated to be 27° (ref 4 ) and 34° 
(ref&). Packing in a crystalline environment planarizes the chain, and in this case the torsional angle is estimated to be 
17° (re^). The quinoid structure of the excited state also planarizes the chain, because in this structure the bridging 
bond has more double bond character, and thus twisting the rings reduces the bond integral and hence increases the 
energy more than in the benzoid structure. The torsional angle in the middle of the distortion is estimated to reduce 
to ~ 9° (ref&). As discussed in section |D] we do not model bond rotations in this work. 

IV. INTERACTING LIMIT 

It is well established that electron-electron interactions enhance the bond alternation in the ground stata 1 ^ and 
generally enhance the size of the lattice distortions for excited states&iS of linear polyenes. The enhancement is 
greater and the electron-lattice relaxation energy is larger for states with covalent character relative to states that 
are entirely ionic in character. Thus, the l 3 ^^ and 2 1 A£ states undergo a greater electron-lattice relaxation than 
the ^B^ u state. These features also occur for the electronic states of light emitting polymers, as we describe in this 
section. First we describe the DMRG algorithm for solving the Pariser-Parr-Pople-Peierls model. 

A. Density Matrix Renormalization Group 

The DMRG method used is basically that of ref 2 , although instead of a three-block method, a four-block method 
is used, with two "middle" blocks being inserted into the middle of the chain each time the chain is grown. 
The procedure is as follows: 

• Initial, exact blocks (spanning the full Hilbert space) are constructed containing three carbon sites (or 7r-orbitals) 
each. These are augmented in the 4-block DMRG program, which creates the biphenyl structure, to construct 
optimized blocks for a phenyl ring. 

• These blocks are further optimized by constructing a four-ring system, with the middle two blocks being obtained 
from a biphenyl calculation with periodic boundary conditions applied. End and middle blocks are reoptimized 
in-situ by performing single-block rotations and truncations. 
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FIG. 5: The ground state energy (relative to the lowest calculated energy) plotted against Hilbert space size for a number of 
different values of m and e. Groups of data points with the same value for e are labelled. 

TABLE II: The vertical and relaxation energies of para-phenylene oligomers (in eV) calculated from the Pariser-Parr-Pople- 
Peierls model. 

State Vertical transition energy Relaxation energy 





N = 4 


N = 8 


N = 4 


N = 8 




3.29 


3.17 


0.58 


0.41 




4.21 


3.96 


0.17 


0.06 


2^+ 


5.52 


5.26 


1.28 


1.14 



• The optimized middle and end phenyl ring basis states are then used in the infinite lattice algorithm to grow 
larger oligomers, with two middle blocks being inserted with each iteration and augmented with the end blocks 
from the last iteration to provide the end blocks for the next iteration. 

• At the target chain size a finite-lattice sweep is performed. 

1. Convergence 

The DMRG convergence parameters used were determined by calculating the ground state energy for a range of 
different values, the results of which are shown in Fig. El There are two truncation parameters in the calculation, m is 
the truncation parameter for the spin and charge sector of a given block with largest number of states. This controls 
the overall truncation of a block. The superblock Hilbert space is controlled by e, the product of the density matrix 
eigenvalues for each block state comprising a superblock state. The size of the Hilbert space is largely independent of 
m and the speed of the calculation is largely determined by the size of the Hilbert space. Since for a given e increasing 
m yields improved results, but has little effect on the speed of the calculation, m was fixed at 50. The results at 
e = 3 x 10 -12 are within approximately 0.01 eV of the lowest energy values, so can be taken to be converged to 
approximately that order of magnitude, and this value was used for all the calculations reported. 

B. Results 

Table HTl lists the vertical and relaxed energies of the 1 3 B^~ U , l x Bi u and 2 1 A+ states for 4 and 8 ring para-phenylene 
oligomers. As in linear polyene s 9 ' 14 , the relaxation energy of the l 1 ^^ state is small, whereas the relaxation energy of 
the 1 3 B^~ U state is large. The experimentally determined relaxation energy of the l 1 !?^ state in the related polymer 
poly(para-phenylene vinylene) has been reported as 0.07 in ref 15 . We may also deduce the relaxation energy in 
poly(para-phenylene) and ladder poly(para-phenylene) from Fig. 3 of refi& by noting that the ratio of the intensities 
of the — 1 to — vibronic peaks in the absorption or emission spectra is S, the effective Huang-Rhys factor. The 
relaxation energy is then fiu x S, where hu) is the characteristic phonon frequency ~ 0.2 eV. Thus, using S — 0.6 for 
ladder poly (para-phenylene), S — 1.2 for poly(para-phenylene) and %uj — 0.2 eV gives relaxation energies of 0.12 eV 
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FIG. 6: The fractional change in transfer integrals of eight-ring (para-phenylene) oligomers from the uniform value, t, in the 
interacting limit. The labels refer to the bonds shown in Fig. [5] Only the upper rung of bonds are shown. 



and 0.24 eV for ladder poly(para-phenylene) and poly (para-phenylene), respectively. The larger relaxation energy for 
poly(para-phenylene) is expected, as the rings are free to rotate, and this result is consistent with a calculated value 
of 0.22 eV reported in ref&. The relaxation energy of the l 1 ^?!^ state in the interacting limit is intermediate between 
the relaxation energy of the state and polaron in the non-interacting limit, as listed in Table [U This illustrates 

the exciton-polaron nature of the l 1 ^?^ state, as further discussed below. 

The relaxation energy of the 2 1 A+ state is also large, but not large enough to cause an energy level reversal of 
the l^-B-j^ and 2 1 A+ states. The difference in relaxation energies between the i 3 B^ u and l 1 ^^ states increases 
the — energy singlet-triplet exchange gap from the vertical gap of ~ 0.6 eV to ~ 0.9 eV, in good agreement 
with experiment*--. We also see that the relaxation energy reduces with chain size, consistent with an increased 
delocalization of the excitations and consequently a diminished effective electron-lattice coupling, and in agreement 
with a wide number of experimental results^. 

Next, we consider the associated geometrical structures. These are plotted in Fig. for the normalized changes in 
transfer integrals and in Fig.0for the staggered, summed bond distortions. As predicted, the ground state alternation 
is enhanced in the interacting limit over the non- interacting limit by 8%. The bond lengths, calculated using Eq. J5J, 
are shown in Fig. [21 

The l^-B-j^ state is now an exciton-polaron**. Its structure is qualitatively similar in both the non- interacting and 
interacting limits, as the soliton-antisoliton confinement due to linear confinement arising from the effective extrinsic 
bond alternation has a rather similar effect to electron-hole attraction. However, as already predicted, the 1 3 B^~ U state 
has a more pronounced distortion because it has some covalent character. Indeed, there is a change of sign in the 
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FIG. 7: The staggered, normalized and summed bond distortions of eight-ring para-phenylene oligomers (as defined in Eq. 
JIZJ) in the interacting limit (filled symbols and solid lines) and non-interacting limit (empty symbols and dashed lines). Note 
that in the non-interacting limit the l 1 -/?-^ and 1 B^ u states are equivalent. 

effective bond alternation. The middle bridging bond becomes a 'short' bond, while the adjacent phenyl-ring become 
'long' bonds. Similarly, the 2}A^ state shows a significant structural distortion, with a change of sign of the bond 
alternation. The lattice distortions of the l 1 !?^, l 3 ^!^ and 2 1 A^ states - as defined by the summed bond distortions 
of Eq. (|17f) and shown in Fig. [7]- are qualitatively similar to those of linear polyenes with extrinsic dimerizarionia. 

The different relaxation energies and geometrical structures of the singlet and triplet B\ u states in the interacting 
limit is obviously related to the different kind of solitons comprising these states, as described in section IIIII In 
particular, the electronic interactions induce a strong coupling of the neutral soliton to the bond-order correlation, 
causing a significant distortion for the triplet state. In contrast, the charged solitons weakly couple to the bond-order 
correlation, and thus the singlet state is more weakly coupled to the lattice. Since the 2}A^ state has an admixture 
of charged and neutral solitons, it also couples more strongly to the lattice than the ^B^ u state. 




V. CONCLUDING REMARKS 



Understanding electron-lattice relaxation, self-trapping and calculating Huang-Rhys factors has important implica- 
tions for predicting electronic processes in conjugated polymers, for example, exciton migration, recombination and 
inter-conversion mechanisms. By solving the Pariser-Parr-Pople-Peierls model of conjugated polymers for poly(para- 
phcnylene) by the DMRG method we have shown that the lattice relaxation of the ^ l B^ u state is quite different from 
that of the 1 3 B^ U and 2 1 A g states. In particular, the ^B^ u state is rather weakly coupled to the lattice and has a 
rather small relaxation energy ~ 0.1 eV. In contrast, the l 3 B^ u and 2 1 A g states are strongly coupled with relaxation 
energies of ~ 0.5 and ~ 1.0 eV, respectively. By analogy to linear polyenes, we argue that this difference can be 
understood by the different kind of solitons present in the l 1 !?^, l 3 Bf u and 2 1 A+ states. The difference in relaxation 
energies of the l x B^ u and 1 3 B^ U states accounts for approximately one-third of the exchange gap in light-emitting 
polymers. 

The results of this calculation present a number of questions. First, are vertical and relaxed states solvated by the 
environment by the same amount ? If not, then the energy gap between the vertical and relaxed states for a polymer 
in the solid state will be different from that predicted here. Second, is the large relaxation energy of the 2}A^ state 
evident experimentally ? Finally, are the different relaxation energies of the \ X B^ U and l 3 ^?^ states relevant to the 
issue of the singlet-triplet exciton fraction in light emitting polymersi2i22i2ii2£ ? 
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